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We calculate the cross section for tjb production at the Tevatron at next-to-leading order in the 
strong coupling and find that more than two millions of rib's are expected per inverse picobarn of 
integrated luminosity. We discuss the decay modes into charmed states and suggest that the decays 
into D* D^*' mesons might be the most promising channels to observe the r/b in Run II. 

PACS numbers: 12.38.Bx,13.25.Gv 



Considering the rich phenomenology of the T 
states, it is quite surprising that spin singlet bb 
states, including the ground state, have not been 
observed yet. 

Fine and hyperfine splittings of the quarko- 
nia spectra were calculated using phenomenological 
models for the heavy-quark potential 0, 0] . Recent 
progress both in lattice and in perturbative QCD 
has allowed the achievement of comparable preci- 
sions 0] • Various approaches have been success- 
fully adopted to describe the charmonium system 
and are believed to be even more reliable for the 
heavier bb states, where relativistic effects are less 
important. Recent determinations lead to a mass 
splitting between the ~£(1S) and the r/h(lS) in the 
40-60 MeV range H S II IS . 

Searches of the rjh have been pursued in various 
experiments. In e~e + collisions, cross sections for 
producing spin singlet states are generally small and 
the signal is rate-limited. This is compensated by a 
clean environment, which allows for searches in the 
inclusive decay modes. Following an original sugges- 
tion by Godfrey and Rosner ||, the CLEO collabo- 
ration has looked for the rjh in the hindered Ml decay 
of the T(35) and found no signal 0. At LEP II, the 
ALEPH collaboration analysed the 77 interactions 
data, basing their a naly sis on the QCD prediction 
of T 77 partial width [lfj : no evidence was found in 
the four- and six-charged-particle decay modes 0] ■ 

The situation is exactly reversed in hadron colli- 
sions, where the production rates can be large, with 
millions of events produced at the Tevatron per in- 
verse picobarn of integrated luminosity, yet the in- 
tense hadronic activity makes any inclusive analysis 
unfeasible. In this case it is necessary to identify 
decay modes that have triggerable signatures and 
allow for full invariant mass reconstruction of the de- 
caying state. Such exclusive modes are believed to 
have very small branching ratios. Braaten et al. ^3 



have suggested that the t% at the Tevatron could be 
observed via its decay into J/ipJ/ip, with the sub- 
sequent leptonic decay of the J/ifi's; experimental 
efforts have started in this direction This sig- 
nature exploits the upgraded abilities of the CDF 
detector of triggering on soft muons. 

The purpose of this note is twofold. First, we 
provide the prediction for the inclusive cross sec- 
tion of rjh production at the Tevatron, based on the 
(NR)QCD calculation at the next-to-leading accu- 
racy in the strong coupling [T"ii . We find that the 
cross section is about four times larger than the 
previous available estimate |l2j . Second, stemming 
from an estimate of the corresponding branching ra- 
tio, we suggest a new analysis based on the decay 
of the rj b -> £>*£>(*) meson pairs. Thanks to sec- 
ondary vertex trigger capabilities of CDF and to the 
high resolution achievable in the invariant mass re- 
construction (~ 20 MeV), this could turn out to be 
the most promising search channel at the Tevatron. 

According to the NRQCD factorization ap- 
proach 0], the inclusive cross section for the rjh in 
pp collisions can be written as: 

a(pp -> Tj h + X) = 

^2 / dx 1 dx 2 f l/p f : j/pa(ij -> rj h ) , (1) 
i,j J 

where 

&(ij^m) = J2 c n( \ o % l°>- ( 2 ) 

n 

The short-distance coefficients , calculable in per- 
turbative QCD, describe the production of a quark- 
antiquark pair state with quantum number n, while 
the (0|OJJ 6 |0) are the non-perturbative matrix ele- 
ments that describe the subsequent hadronization of 
the bb pair into the physical rjh state. These matrix 
elements can be expanded in powers of v 2 ~ 0.1, the 
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FIG. 1: Representative Feynman diagrams for r/b hadro- 
production at LO (a), and virtual (b) and real (c) con- 
tributions at NLO. 

relative velocity of heavy quarks in the bound state, 
so that, to a given accuracy, only a few terms need 
to be included in the sum over n. 

The case of gb production is particularly simple. 
Compared to the leading contribution, the 1 So~\bb), 
the color octect terms, ^Sf 1 , 3 Sf 1 , 1 Pf 1 ), are all 
suppressed by u 4 , while the corresponding short dis- 
tance coefficients start at least at a| as the sin- 
glet contribution j^. Hence, given that the non- 
perturbative matrix elements for the singlet produc- 
tion are extracted from T decays, there are no un- 
known parameters (up to corrections of 0(v )) en- 
tering the estimate for i% production. 

The leading-order cross section for gg — > rjb, 
Fig. da), is 

^99 - m) = - ^){o\o?Cs )\o). (3) 

Som^s s 

Next-to-leading order corrections in the strong cou- 
pling have been calculated in Rcfs. Hi [3. These 
include virtual corrections to the 2 — ► 1 process, 
Fig. ^b), and real 2^2 processes, such as gg — > 
rjbg,gq — > rjbq, Fig. ^c). The result for the total 
cross section in pp collisions at 1.96 TeV of center- 
of-mass energy is 

a(pp -> Vb + X) = 2.5 ±0.3 ^b, (4) 

where we have adopted CTEQ5M1 parton densities, 
the corresponding two-loop evolution for as(pn) 
and mb = 4.75 GeV. The quoted uncertainty has 
been estimated by summing (in quadrature) the er- 
rors coming from various sources. The first is as- 
sociated to the choice of the renormalization and 
factorization scales. Varying them independently in 
the range m;, < /ip < Arrib gives an uncertainty 
of 10%. For the non-perturbative matrix element we 
have used the determination from T leptonic decay, 
(T|C 1 ( 3 S' 1 )|T} = 3.5 ± 0.3 GeV 3 , which can be re- 
lated to the rjb production matrix element by using 
spin-symmetry and vacuum saturation approxima- 
tion, (0\O? (%)|0) = (T|Oi( s 5i)|T>, up to 0(« 4 ) 
corrections. Other determinations coming from po- 
tential models Il7j and lattice calculations [18| fall 



in the range of the quoted error. Finally, we also 
included the effect of an uncertainty on the bottom 
(pole) mass of ±50 MeV on the cross section. The 
strong correlation to the non-perturbative matrix el- 
ement extraction has been exploited to reduce the 
uncertainty from this source. The result, quoted 
in Eq. 0}, accounts for the "direct" contribution 
and does not include feed-downs from higher-mass 
states, such as the hb- 

For the experimental analysis it is important to 
know the distribution of r)b at small p T values. 
This cannot be described accurately by a fixed-order 
calculation. In this region of the phase space it 
is necessary to resum higher-order corrections in- 
volving soft-gluon radiation. To this aim, we use 
PYTHIA [lj|, matched with the exact matrix el- 
ements for ij — > ri p k describing the high-p T tail, 
as outlined in Ref. [20(. This procedure, which has 
been shown to work well for the analogous process 
gg —> has the additional virtue that it can be 
directly used for simulation in experiments. The 
results are shown in Fig. [21 where the differential 
distribution in p T for the r)b is shown (upper curve). 
Hadronization and initial kx effects are not included. 
The normalization of the inclusive r\b distributions is 
obtained from the NLO result, Eq. 0}. 

As already stated above, the cross section for the 
rib in the central region is about four times larger 
than was estimated in Ref. ^| by rescaling the T 
cross section at high transverse momentum. This 
procedure is expected to provide an underestimate, 
since it assumes that the p T spectra for T and rjb 
have a similar shape at small p T . In fact T color- 
singlet production proceeds at LO through gg — » Tg 
and therefore vanishes at p T ~ 0, where the largest 
part of the rjbS are produced. 

With such a large number of events expected, it is 
interesting to consider in detail the rare exclusive de- 
cays that might give triggerable signatures. It is easy 
to show that direct decays into photons or lepton 
pairs give either too small branching ratios or very 
difficult experimental signatures (such as rip — » 77 
whose branching ratio is O(10 -5 )). In Ref. [12|, the 
branching ratio Br(r/b — ► J/ipJ/ip) was estimated 
through scaling from the analogous r\ c — ► (jxfi, finding 
a value compatible with 7 x 10~ 4±1 . This is proba- 
bly an overestimate. As a simple upper bound, let 
us consider the inclusive decay rate of the 77b into 
four-charm states 

r(% -> J/i/i J /ii)) < r{ Vb -> cccc). (5) 

The inclusive rate can be calculated at leading order 
by considering the four Feynman diagrams, such as 




pp -> 77 b + X @ 1.96 TeV 




FIG. 2: Differential cross sections for 775 production at 
the Tevatron (pp at 1.96 TeV). The upper curve is nor- 
malized to the NLO total rate and describes the p T dis- 
tribution for the r/t- The lower curve is the correspond- 
ing distribution for the D mesons coming from the rjb 
decays, after requiring that they both have \rj(D)\ < 1 
(no branching ratio included). After the acceptance cut, 
the rate drops to 15% of the total cross section. 

the one shown in Fig.JSfb). The result is 

Br(7/- 6 -> cccc) = 1.8±g;l x 1(T 5 . (6) 

where m c = 1.45 GeV, nib — 4.75 GeV, a>s(2mb) = 
0.182 and the NLO expression for the inclusive to- 
tal width |3 H3 have been used. The amplitudes 
have been calculated analytically while the integra- 
tion over the phase space has been performed numer- 
ically. The uncertainty has been estimated by vary- 
ing the renormalization scale between nib < /xr < 
477ifc and the masses in the ±50 MeV range. The 
four-charm branching ratio is very sensitive to the 
value of the charm mass, which dominates its uncer- 
tainty. The above result shows that the inclusive 
rate is already smaller than the lower bound ob- 
tained in Ref. |12J; further suppression is expected 
mainly because many other decay modes to charmed 
mesons (or other charmonium states) should con- 
tribute to the saturation of the inclusive rate. 

Given the result in Eq. 0, a comment on the re- 
liability of the estimate based on the scaling from 
Br(77 c — * (jxf) is in order. To this aim we recall 
that the decay of a scalar QQ meson into two vec- 
tor states is suppressed in perturbative QCD [2^ . 
A non-trivial check of this selection rule is that 
the branching ratio Br(r7;, — > J/tpJ/tp) is exactly 
zero when calculated at LO in the NRQCD dou- 
ble expansion in as and v 2 . For the same reason 
one would expect the rate of r\ c — > (fxp to be sup- 
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FIG. 3: Representative Feynman diagrams for the rjb 
inclusive decay into (a) two-charm (two diagrams) and 
(b) four-charm (four diagrams) states. 

pressed, in contradiction with the measured value 
of about 1%. This entails that some other (non- 
perturbative) mechanism is responsible for this de- 
cay process |2j|. Rescaling by (m c /m,b) 4 the branch- 
ing ratio of n c — > 4>4> to obtain the branching ratio 
of r]b — * J/ipJ/ip amounts to rescaling by the same 
factor also the effect of non-perturbative or higher- 
order contributions that are likely to be crucial in 
determining the rj c decay, but less and less impor- 
tant as we pass from the r] c to the 775 system. 

From the phenomcnological point of view, Eq. © 
implies Br(?7fc — > /i + /i~/i + /i~) < 6 x 10~ 8 , which 
makes a search in this channel very hard. As an 
alternative, we propose to consider the decays into 
charmed mesons, % —*■ £)*£)(*) |2sj . 

A computation of the exclusive decay rates is 
missing and does not seem to be feasible within 
the framework of quark models or QCD sum rules. 
Probably, there is room to face such a problem using 
lattice techniques. Here we assume that the exclu- 
sive decays into D*D^ dominate the inclusive rate 
into charm: 

T{n b ^ D*D^) <T(r lb ^ cc + X). (7) 

We find that the largest contribution to the decay of 
the r/b into charmed states is given, Fig. 1ST a), by 

Br(Tfe->ccg) = 1.518:2%, (8) 

where we used the same input parameters and es- 
timated the uncertainties as in the decay into the 
four-charm states. It can now be argued that not 
much suppression is expected in passing from the 
inclusive to the exclusive decays. For example, 
decay rates into charmonium states should be of 
the same size T -► J/ip + X, i.e., O(10" 3 ) 0. 
Emission of extra pions should also be considered, 
but some of these contributions would be automat- 
ically included in the experimental analysis, e.g., 
rjb — » DD* — » DDir, as non-resonant diagrams. 
However, to be conservative, we consider the range 
10~ 3 < Br(77 fc -> D*7J(*)) < 10~ 2 in the following 
phenomcnological analysis. 
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There are other decay modes leading to a two- 
charm final state that have not been included in our 
estimate of the branching ratio. One is the decay of 
rjb — > g* g* — > cc through a (box) loop, which pro- 
ceeds at order a A s and is further suppressed by loop 
factors. Another is the decay rjb — g* — > cc, via 
a 3 5'[ 8 ' state, which is only, a 2 s but it is suppressed 
by the color octet matrix element (7jb\Os( 3 Si)\rjb). 
Assuming a scaling 0(v ) for the non-perturbative 
matrix element, this process gives a non- negligible 
contribution to the branching ratio, about 5 x 10~ 3 . 
However, this result is affected by a large uncer- 
tainty and could be much smaller as suggested by 
studies I25L l26j about the size of color octect matrix 
elements in T decays. 

Results for the p T distribution of the D mesons 
from the rjb decay are shown in Fig. (lower curve). 
No branching ratio is included, the difference in 
rate coming only from the requirement that both 
D mesons be central, \rj(D)\ < 1. As expected, the 
p T distribution peaks just below M Vb /2. The effi- 
ciency for the geometrical acceptance of the detector 
is found to be about 15%. By adding the require- 
ment that at least one D meson has p T > 5 GeV, one 
is left with only 4% of the total number of events pro- 
duced. The above efficiency can be folded with our 
estimate of the branching ratio leading to 10 4 — 10 5 
£)*£>(*) triggerable events expected from the rjb de- 
cay in 100 pb _1 of integrated luminosity at the Teva- 
tron. The final number of reconstructed events will 
depend on the decay modes of the D and D* mesons 
and on the associated experimental efficiencies. We 
leave this to more detailed experimental studies and 
only add a few comments. First we note that, ac- 
cording to the arguments outlined above, perturba- 
tive QCD predicts the rib — ► D*D* decay to have 
a smaller rate with respect to rjb — > D*D. In this 
case, it is reasonable to expect that different charge 
assignments, such as D*° D° , D*+ D~ , D*~ D+ , will 
occur with the same probability of ^ . Finally, we re- 
call that the cleanest signatures have small branch- 
ing ratios, Br(D° K~tt + ) = 3.90% and Br(D + -> 

K n + ) = 2.77% [13], leading to a factor of about 
10~ 3 drop in the rate if both D mesons are required 
to decay through these channels. We foresee that 
sizeable improvements could be achieved by requir- 
ing that just one of two D mesons decays through a 
very clean signature, providing an efficient trigger. 

To summarize, we have presented the NLO QCD, 
prediction for rjb production at the Tevatron, includ- 
ing a resummed result for the p T spectrum obtained 
with a dedicated implementation in PYTHIA. The 
production rate is large, of the order of a few jib, 
and allows for the search of the rjb through rare ex- 



clusive decays. We argued that the branching ratio 
into J/ip J/ip is probably too small and we suggested 
to look for the rjb through its decay into D*D^ 
mesons. Our results indicate that the rjb could be 
eventually observed during Run II at the Tevatron. 
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